Modeling rhythmic patterns in the hippocampus 
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We investigate different dynamical regimes of neuronal network in the CA3 area of the hip- 
pocampus. The proposed neuronal circuit includes two fast- and two slowly-spiking cells which are 
interconnected by means of dynamical synapses. On the individual level, each neuron is modeled by 
FitzHugh-Nagumo equations. Three basic rhythmic patterns are observed: gamma-rhythm in which 
the fast neurons are uniformly spiking, theta-rhythm in which the individual spikes are separated by 
quiet epochs, and theta/gamma rhythm with repeated patches of spikes. We analyze the influence 
of asymmetry of synaptic strengths on the synchronization in the network and demonstrate that 
strong asymmetry reduces the variety of available dynamical states. The model network exhibits 
multistability; this results in occurrence of hysteresis in dependence on the conductances of indi- 
vidual connections. We show that switching between different rhythmic patterns in the network 
depends on the degree of synchronization between the slow cells. 

PACS numbers: 87.18.Sn, 87.19.11 
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I. INTRODUCTION 

Rhythmic behavior is a basic property of biological sys- 
tems and plays an important role in a variety of physi- 
ological processes In particular, hippocampal neu- 
rons of the human brain are able to generate rhythmic 
oscillations in several frequency ranges; among these, 
prominent roles belong to fast (gamma-rhythm, 30^-65 
Hz) and slow (theta-rhythm, 3-^10 Hz) signals and to 
mixed regimes in which fast and slow oscillations alter- 
nate. Theta-rhythm has been demonstrated to be re- 
sponsible for coding of spatial information [2j, synaptic 
modification of intrahippocampal pathways, as well as 
for assembling and segregation of neuronal groups Q. 
In its turn, oscillatory activity in the gamma-frequency 
band is involved in information transmission and stor- 
age [|||. Isolated parts (so-called CAl, CA3, DG etc.) 
of hippocampal formation where these regimes are ob- 
served, have been subjected to numerous experiments 
foT-flfjj. In particular, it has been established that the 
neuronal network in the CA3 region of the hippocampus 
can exhibit three qualitatively different oscillatory states 
(theta-, gamma- and the mixed "theta/gamma") and is 
able to dynamically switch between them Q ■ Transitions 
between different regimes in brain neuronal networks can 
be provoked not only by variation of parameters, but by 
variation of initial conditions as well [TT14T3 ; the latter 
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may occur due to changes in the baseline potential [13j 
or to perturbations of initial ionic concentrations (l5j . 

Cells, which constitute a hippocampal network, differ 
in their morphologi cal, electrophysiological and neuro- 
chemical properties [lq . Synaptic interactions between 
several types of cells: the pyramidal ones, basket, sep- 
tum and oriens-lacunosum-moleculare (OLM) cells are 
responsible for the generation of different rhythms 10]. 

Besides experimental studies, transitions between var- 
ious rhythms and multistability phenomena in isolated 
cells as well as in neuronal networks have been subjected 
to theoretical and numerical modeling @, fl7i - [20j . Typ- 
ically, such models involve various cell types as well as 
large numbers of ionic currents and compartments. In 
the detailed models (l9l l2lll as well as in simpler bio- 
physical approximations [17( the cells, responsible for the 
generation of rhythmic patterns, are typically described 
by equations of the Hodgkin-Huxley type. 

Role of intercellular synchronization in the transitional 
processes is an object of ongoing discussions. In particu- 
lar, it has been shown [22[ that a switching between two 
states of the system can forecast desynchronization in 
a bursting activity of the hippocampus. The strength of 
intercellular coupling governed by synaptic conductances 
is apparently a crucial factor in this context @, HH ■ For 
example, studies of detailed thalamocortical model [20j 
have shown that presence or absence of bistability in the 
neuronal circuit can depend on parameters of the electri- 
cal coupling between neurons as well as on the synaptic 
conductance. 

High degree of complexity, obligatory in realistic mod- 
els, hampers revealing of principal mechanisms and con- 
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trol parameters of transitions between different types of 
oscillations. Therefore, along with accurate reproduction 
of experimental observations, one of the main aims of the 
modeling remains a consideration of general dynamics in 
the appropriate class of circuits, as well as an understand- 
ing of the underlying dynamical mechanisms which are 
responsible for different rhythmic patterns. This raises 
a question about a minimal module of the hippocampal 
network: a small set of cells capable of reproducing all ba- 
sic typs of experimentally observable rhythms [17]. Em- 
bedded into intricate networks of neurons, such modules 
can play a role of pacemakers which ensure persistence 
of optimal oscillatory states in the entire ensemble. 

Below we discuss a minimalistic dynamical model of 
the hippocampal circuit in the area CA3, which is able 
to produce several rhythmic patterns. For this purpose 
we simplify the basic module of the circuit in such a 
way that maintenance of the essential types of dynami- 
cal regimes would be controlled by a concise number of 
principal parameters. Since we do not aim at detailed 
reproduction of intrinsic cellular dynamics but restrict 
ourselves to rhythmic aspects, we reduce the description 
of each element to the classical simplified model of neu- 
ronal oscillations: the FitzHugh-Nagumo (FHN) equa- 
tions. Simplification of elements in this case does not 
imply trivialization of dynamics: ensembles of the FHN 
oscillators are known to exhibit different kinds of col- 
lective behavior which, depending on the topology and 
strength of interaction, range from subthreshold oscilla- 
tions via mixed- mode intermittency to tonic spiking [2~ij . 
Interplay between connectivity and coupling strength in 
networks of heterogeneous FHN oscillators can induce 
or destroy synchrony [25| . Working with FHN neurons 
renders the well-understood dynamics on the level of in- 
dividual elements and allows us to concentrate on the 
governing role of interactions between them. 

In Sect. [IT] we introduce the model equations and dis- 
cuss the choice of appropriate parameter values. Sect. IIIII 
begins with characterization of the principal rhythmic 
patterns; in particular, the role of the phase shift be- 
tween the OLM cells is discussed. Further we proceed to 
the description of multistability and transitions between 
the patterns. We demonstrate that the model repro- 
duces the basic rhythms observed experimentally, as well 
as hysteresis between them. Finally, we show how the 
symmetry or asymmetry in the arrangement of synap- 
tic connections influences the existence and properties of 
dynamical states. 



II. MODEL EQUATIONS 

In a recent paper [7( , a kind of minimalistic network for 
studies of different regimes in the hippocampal area CA3 
has been proposed. This network includes main types 
of hippocampal cells and consists of five elements: two 
fast-oscillating cells (so-called basket cells) , two slow cells 
(OLM - oriens/lacunosum-moleculare associated - cells) 



and one two-compartmental pyramidal cell. The cells are 
synaptically connected in all-to-all topology, with one ex- 
ception: there is no direct connection between the slowly 
spiking cells. Accordingly, the pyramidal cell activates 
the rest of the cells which, in turn, inhibit each other 
and the pyramidal cell. The cells are described within 
the framework of the Hodgkin-Huxley formalism. After 
taking into account all relevant currents as well as kinet- 
ics of synaptic variables, the resulting system constitutes 
a set of 41 ordinary differential equations. Notably, the 
choice of coupling coefficients employed in Q ensured 
a certain degree of asymmetry in the network: in par- 
ticular, each basket cell was coupled to both OLM-cells 
through connections with different conductivity. 

Numerical simulations of that model have confirmed 
that transitions between different regimes depend on the 
strength of coupling between the cells. Typically, a 
switching between different rhythms is preceded by on- 
set of a certain phase shift between two slow-spiking cells. 
Within the model of @, this shift can be ensured e.g. by 
a judicious choice of appropriate initial conditions. 

In the present work we further reduce the complex- 
ity of the model: we decrease the number of cells and 
replace the Hodgkin-Huxley setup by the FitzHugh- 
Nagumo equations. The resulting network, sketched in 
FiglTJ includes one-compartmental pyramidal cell P, one 
basket cell B and two OLM-cells L\ and L^. All cells 
are mutually synaptically connected; the only exception 
are two slow OLM-cells which do not communicate di- 
rectly. Synaptic connections are governed by the synaptic 
conductances Gji with j, i = L\, L 2 , P, B (and Gjj=0). 
Leaving in the network only one basket cell B, we take 
into account, however, inequality of the cross- and direct- 
connections as a potential source of a phase shift between 
the slowly spiking cells. For this purpose, it is sufficient 
to allow for asymmetric synaptic inputs from the basket 
cell to L\ and L^. in general, GpLt ^ Gbl 2 - The rest 
of connections is set symmetrically: Gpl 1 =Gpl 2 = Gpl, 

LiP =< ^L 2 P =, ^LP, ( ~'L 1 B=t'L 2 B= ( ~'LB- 

Within this simplified description, it appears reason- 
able to keep only one activating input which acts upon 
the fast-oscillating cell P. The latter excites all other 
cells, which inhibit each other as well as P itself. 

We use below the form of equations in which all vari- 
ables, and, hence, all parameters are measured in dimen- 
sionless units. Each cell of the network obeys the stan- 
dard set of coupled FitzHugh-Nagumo equations: 

i 



where Vi is the membrane potential of the i-th cell (recall 
that i = L%, L12, P, B), and ut is the respective membrane 
variable. The values of parameters a and b are the same 
for all cells. The Kronecker-delta Si p ensures that only 
the pyramidal cell is externally excited by the current 
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FIG. 1: (Color Online). Sketch of the cellular network. L\^'- 
slow OLM cells; P: fast pyramidal cell; B: basket cell. Ar- 
rows: directions of currents. Filled circles: synapses. Pluses: 
excitatory synapses; minuses: inhibitory synapses. Synaptic 
connections are characterized by conductances Gji. External 
input Iext affects only the pyramidal cell P. 



The synaptic input /|yn in Eqs (JlJ - the current from 
cell j to cell i - is defined as 

/CrO =G--s-(E ■ -v) 

and is governed by kinetic equation for the synaptic vari- 
able Sji(t) H3: 

^gi = ±(l + tanh IL) (1 - Sji ) - B Sji , (2) 
at l \ v vs i / 

For simplicity, we assume that all synapses obey quan- 
titatively the same kinetics. This is achieved by using 
the same values of A, B and v vs i in © for all Sji. Hence, 
evolution of Sji depends (via Vj) only on j, and the num- 
ber of independent synaptic variables equals the number 
of cells in the network. Altogether, the dynamical system 
includes 12 ordinary differential equations: 8 equations 
for the FitzHugh-Nagumo variables along with 4 equa- 
tions for the synaptic variables. 

Originally, the FitzHugh-Nagumo equations were writ- 
ten as a crude simplistic model which mimics the charac- 
teristic features of the Hodgkin-Huxley dynamics. Hence, 
their parameters do not allow for direct biological inter- 
pretation and serve the mere purpose of supporting the 
proper dynamical regimes. To illustrate possible kinds of 
dynamics in the considered network, we fix the parame- 
ters of the FitzHugh-Nagumo equations ((T|) at the values 
which ensure that in absence of external currents every 
isolated cell is at rest; if perturbed by sufficient exter- 
nal input, it displays full-scale oscillations. The values 
a=0.5, 6=0.8 which we use for our calculations, are close 
to those employed by FitzHugh |26|. 



The parameter Si in the FitzHugh-Nagumo equations 
(U} determines the timescale of oscillations in the isolated 
i-th cell: in the leading order, period of relaxation oscil- 
lations in a cell is inversely proportional to £j. The di- 
mensionless number £j can be interpreted as ratio of two 
dimensional timescales of the neuron: that of the fastest 
relaxation to that of the slowest one. In the case of the 
pyramidal and the basket cells these are the relaxation 
times of, respectively, potassium and sodium currents; 
the values of £j for these two cells can be viewed as (ap- 
proximately) the same. In contrast, in the OLM-cells 
the slowest processes involve relaxation of specific so- 
called hyperpolarization-activated currents and 
their values of er, should be chosen substantially lower. 
Accordingly, we set £=0.3 for the fast basket and pyra- 
midal cells, whereas for slow L-cells the value e=0.04 is 
adopted. 

Equations ^ for synaptic variables are phenomeno- 
logical as well; the values of parameters A, B and v vs i 
are commonly obtained by fitting the experimentally 
measured dependencies. By setting A=l, B=0.3 and 
v vs i=0.l we model the typical situation in which opening 
of channels occurs faster - but not much faster - than 
their subsequent closure. 

The value of E in the expression for synaptic current 
depends on the kind of connection. For all excitatory 
synapses we put E cx = 0, whereas for inhibitory synapses 
Ein=~5 is set. The latter value ensures that the difference 
E m — Vi , and thereby the corresponding synaptic current 

(ji) 

Jsyn stay negative at all times, and the connection is 
indeed inhibitory. 

Concerning synaptic conductances Gji, we choose 
the values which, dimensionahzed in units of 
0.01O/iTO _1 m -2 , have the order of magnitude of 
those, reported in Q- In a sense, conductances prede- 
termine the relative importance of individual cells in 
the ensemble. Since external current acts only upon the 
pyramidal cell, synaptic connections from this cell to the 
rest of the module should be strong enough in order to 
activate the otherwise silent basket and OLM-cells; this 
is ensured by assigning to them relatively high values 
Gpl = 0.7 and Gpb = 0.57. A lower conductivity is 
assigned to the backward connection from the basket 
to the pyramidal cell: Gsp=0.1. Further, we introduce 
asymmetry between the left and right halves of the 
network by fixing different values for connections from 
the basket to the OLM-cells: GbLi=0.06 in contrast to 
G BL2 =0.03. 

The slow OLM-cells are driven by the pyramidal cell; 
intensity of their inhibitory feedback to the driving ele- 
ment depends on conductances Glp and (mediated by 
the basket cell) on conductances Glb- If both Glb and 
Glp are set to zero, the feedback is absent, and the OLM- 
cells play the passive role. Each of them oscillates period- 
ically: it is phase-locked to the oscillations of the driving 
fast subsystem formed by the pyramidal cell and the bas- 
ket cell. However, due to asymmetry, the slow cells can 
be locked (and under the quoted parameter values are 
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indeed locked) to the fast subsystem in different locking 
ratios. Asymmetry of rhythmic patterns persists for suf- 
ficiently low values of Glb and Glp- this is visualized 
in the lower panel of Fig J5] where the slow cells oscillate 
in the ratio 2:3. Increase of conductances Glb and/or 
Glp intensifies the feedback; in spite of the absence of 
the direct connection, the OLM-cells interact through the 
mediation of the fast cells, and their periods get adjusted. 

In order to concentrate on the role of connection from 
the OLM-cells to the pyramidal cell, we use throughout 
this work the conductance Glp as a controlling param- 
eter and study transitions which occur in the course of 
its variation. As for the connection from the OLM to 
the basket cell, it should be weak enough in order not 
to damp the module dynamics, hence we fix for our cal- 
culations (with the only exception of the plot in FigJ^]) 
the value Glb = 0.01. This value is sufficiently high 
to ensure that the feedback via the basket cell synchro- 
nizes the slow cells in the frequency ratio 1:1, even in the 
absence of direct feedback through the LP-connections. 

The steady solution of equations corresponds to the 
quiescent non-spiking state. Intensity of the constant 
external input I cxt (recall: it affects only the P-cell) 
should ensure destabilization of the equilibrium and sus- 
tainment of oscillations in the system. Under Glp=0.037 
and the aforementioned fixed values of other parameters, 
the equilibrium is unstable for 0.02 < I cxt < 1.1. In the 
endpoints of this interval, subcritical Hopf bifurcations 
take place. Location of the endpoints is almost insensi- 
tive to variation of Glp- Accordingly, for our numerical 
simulations we take the value of input from this range: 
4xt = 0.43. 

Numerical solutions have been obtained using the soft- 
ware packages XPPAUT and MATCONT. 



III. RESULTS 

In order to characterize different attractors of the sys- 
tem and to understand possible mechanisms of switching 
between the regimes, we have studied the response of 
the model to variation of the strength of LP connection. 
Starting from the situation when this connection is ab- 
sent, we describe below the states which are observed in 
the course of increase of Glp- 

Among the solutions, we single out the three main 
qualitative types of oscillations. The difference between 
them lies in the behavior of the voltage variables of fast 
(pyramidal or basket) neurons: in the gamma-rhythm 
the voltage exhibits fast regular oscillations, whose am- 
plitude is weakly modulated due to interaction with slow 
counterparts. In the theta/gamma rhythm certain spikes 
in the pattern "fall out" ; the plot shows regular patches 
of spikes interrupted by nearly quiescent epochs. Typ- 
ically, theta/gamma oscillations are periodic. Finally, 
the theta pattern consists of repeated solitary spikes on 
the nearly quiescent background. For a slow OLM-cell 
the difference in the temporal dynamics between three 



rhythms is not so pronounced; here, transitions manifest 
themselves in the onset or the disappearance of a phase 
shift between two such cells. 



A. Gamma-rhythm 

If the value of conductance Glp is set to zero, the 
OLM-cells exert no inhibiting action upon the pyramidal 
cell; this corresponds to removal of the upper connecting 
arcs in the scheme of FigJT] In this situation, the sys- 
tem performs periodic oscillations of the gamma type. 
Consecutive spikes in the fast cells occur at close time 
intervals, but the exact repetition (in other words, clo- 
sure of the limit cycle in the phase space of the system) 
requires several spikes. In particular, smallness of Glb 
can result in asynchrony between the slow cells (Fig. [2]). 
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FIG. 2: (Color Online). Gamma-rhythm at Glp=0.032, 
Glb— 0-001. Top: pyramidal cell; periodic pattern which 
consists of 9 spikes. Bottom row: OLM cells; solid curve: 
Li; dashed curve: 1/2- OLM cells are locked in the frequency 
ratio 2:3. 

By continuity, rhythmic pattern of the gamma type 
persists at sufficiently low non-zero values of Glp as well. 
An example is shown in Fig|3l At Glp = 0.035 the slow 
cells oscillate in the ratio 1 : 1 (bottom panel), but their 
maxima do not occur simultaneously. In this context it is 
natural to speak about the phase shift. We introduce the 
phase phcnomcnologically: for each cell the phase incre- 
ment of 27T is assigned to every interval between two con- 
secutive positive maxima of voltage; between the max- 
ima, the phase is linearly interpolated. The phase shift 
between oscillations shown in the bottom panel of Fig. [3] 
is close to 2-7r/3. 

B. Theta/gamma rhythm 

In a broad range of values Glp above 0.0362 an ap- 
propriate choice of initial condition results in a differ- 
ent rhythmic pattern: the system exhibits theta/gamma 
oscillations. The bottom panel of Fig. 0] indicates that 
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FIG. 3: (Color Online). Gamma-rhythm at G L p=0.035, 
Glb— 0-01. Top: pyramidal cell. Bottom row: OLM cells; 
solid curve: Li; dashed curve: 1/2- OLM-cells are locked in 
1:1 ratio; note phase shift between them. 



in the theta/gamma state both slow cells oscillate syn- 
chronously; there is no phase shift between them. 
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FIG. 4: (Color Online). Theta/gamma-rhythm close to its 
onset; Gip=0.05. Top: pyramidal cell. Bottom: OLM cells; 
solid curve: Ly, dashed curve: Li- 



C. Theta-rhythm 

At sufficiently high values of Glp theta/gamma 
rhythm is replaced by theta-oscillations. Switching be- 
tween the two regimes is visualized in Fig. [5] the motion 
starts as a typical theta/gamma oscillation, with patches 
of spikes of the fast cell; in the course of time evolution, 
at around t w 200, it abruptly changes the pattern and 
acquires the characteristic shape of theta rhythm with 
equidistant solitary spikes of the fast cell. Remarkably, 
the slow cells, synchronous on the initial stage, develop 
during the switch the phase shift of 7r: in the theta- 
rhythm the OLM cells oscillate in antiphase. 
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FIG. 5: (Color Online). Breakdown of theta/gamma rhythm 
and onset of theta rhythm, Glp~2.274. Top: pyramidal cell. 
Bottom: OLM cells; solid curve: L\ \ dashed curve: Li- Note 
in-phase oscillations before the switch and antiphase immedi- 
ately after it. 



D. Hysteresis 

Over large intervals of values of Glp, different rhyth- 
mic patterns coexist as attractors of Eqs (|H2p : at 
low values of Glp, depending on the initial values 
of the variables, the system can exhibit gamma- as 
well as theta/gamma-oscillations, whereas at moderate 
and high values of Glp there is a hysteresis between 
theta/gamma- and pure theta- patterns. 

For slow OLM cells, neither a transition between 
gamma- and theta/gamma-states, nor the further growth 
of Glp significantly affect the period of their oscillations. 
In contrast, periodicity in spiking patterns of the fast 
cells P and B changes quite noticeably. A convenient 
characteristics is delivered by the average duration of in- 
terspike interval (ISI): the mean distance between the 
positive maxima of voltage in a cell. In Fig. [6j we plot 
this characteristics for membrane potential oscillations in 
the pyramidal cell. The value of ISI increases from 14.73 
in the gamma state at Glp = to 35.94 in the theta 
state at Glp = 3. Similar situation is observed for the 
basket cell. 

A closer look shows that gamma- and theta- branches 
are connected; in fact, they belong to the same con- 
tinuous family of solutions. A transformation from the 
gamma- to theta-rhythm occurs through the intermedi- 
ate regime which is observed within the small parame- 
ter range 0.0724 < G L p < 0.0734. Main stages of this 
process, in terms of the voltage variable of the pyrami- 
dal cell, are presented in Fig. [7l The transformation 
is preceded by a gradual deformation of the gamma- 
pattern: of the three initially nearly equal spikes (cf. 
Fig. [3]), two spikes get diminished (Fig. (7£a,b)). As a 
result, the periodic gamma-state acquires the character- 
istic shape of mixed-mode oscillations [28| with alternat- 
ing small subthreshold - and large-scale spiking maxima. 
At Glp — 0.0726 these oscillations form the complicated 
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FIG. 6: (Color Online). Duration of interspike intervals in 
dependence on Glp for different rhythmic patterns. Solid 
line: gamma-rhythm. Dotted line: theta/gamma pattern. 
Dashed line: theta rhythm. Inlet: range of low values of the 
conductance Glp- 



(yet periodic, with period 541.1 !) temporal pattern, vi- 
sualized in Fig. [TJc). A small increase of Glp leads to 
partial flattening of subthreshold epochs (Fig. |7£d)); the 
next small increment results in the regular alteration of 
a spike and a single subthreshold oscillation (Fig. E^e)). 
In fact, this state is already very close to the theta- 
pattern; further increase of Glp brings only quantita- 
tive changes: the subthreshold oscillations are gradually 
flattened, eliminated, and the theta state is established. 
On the phase portrait (cf. Fig. |8|), the epochs of sub- 
threshold oscillations are represented by minor loops of 
the trajectory. In the course of flattening, these loops 
get smaller, turn into cusps and disappear, so that only 
the large loops (corresponding to solitary spikes of the 
theta-rhythm) persist. 

Remarkably, the described transformation changes the 
phase difference between the slow cells: as already men- 
tioned, in the gamma-state the difference is close to 2tt/3, 
and the proximity to this value holds until the onset of 
mixed-mode oscillations in Fig.[71[c). As soon as the sim- 
ple temporal pattern is regained (Fig.Efe)) the phase dif- 
ference acquires the value close to 7r, and the antiphase 
oscillations of slow cells remain the hallmark of the theta- 
state everywhere in the range of its existence. 

A similar process of deformation occurs on the 
other branch of solutions, which corresponds to the 
theta/gamma state. In this state, patches of full-scale 
spikes are separated by subthreshold oscillations (cf. 
FigHJ; the latter correspond to minor loops on the pro- 
jections of respective phase portraits. In the course of 
increase of conductance Glp, the shape of the attracting 
orbit is subjected to continuous deformation (similarly to 
that shown in Fig. |8]for theta-rhythm): the minor loop 
gradually gets smaller, shrinks into a cusp and disappears 
completely. This process leads to gradual decrease and 
elimination of one of the maxima in the temporal pattern 
of the oscillations. 

As seen in Fig. [51 the branch of theta/gamma os- 
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FIG. 7: (Color Online). Transformation of gamma- into 
the theta-rhythm for the membrane voltage of the pyrami- 
dal cell, (a) G L p=0.07; (b) G L p=0.0724; (c) Glp=0.0726; 
(d) Glp=0.0732; (e) Glp=0.0734; (f) Glp=0.1000. 



cillations is isolated in the parameter space. A closer 
look shows that its birth at low values of Glp be- 
gins with a very short segment of mixed-mode oscilla- 
tions; they are present close to Glp— 0.0362, but al- 
ready at Glp=0.0363 the characteristic shape of periodic 
theta/gamma oscillations is established. On the other 
end of the branch, at Glp— 2.274, the limit cycle which 
corresponds to the theta/gamma rhythmic pattern, dis- 
appears in the saddle-node bifurcation. 



E. Attraction basins of coexisting states 

In further numerical simulations we studied the sensi- 
tivity of the system response with respect to a variation 
of initial conditions. At low (below 0.036) values of con- 
ductance Glp the oscillatory state of the gamma-type 
was reached for all tested initial conditions. On the other 
end, at sufficiently large (above 2.28) values of Glp the 
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FIG. 8: (Color Online). Evolution of attractor shape in the 
phase space. Theta-rhythm. Projection on the variables of 
the pyramidal cell. Minor loops, which correspond to sub- 
threshold oscillations, shrink, turn into cusps and disappear. 
Solid line: Glp=0.09, dashed line: Glp=0.2, dotted line: 
Glp=0A. 

system appears to possess the unique attractor as well: 
now this is the limit cycle which corresponds to the theta 
regime. In between, the system is bistable. 

For a case study we take the value Glp = 0.8, well 
inside the range of this parameter where the distinct os- 
cillations of the theta type coexist with the theta/gamma 
rhythmic pattern. 

The system has been scanned on a grid of initial 
conditions for variables Vi and Ui from -2 to 2: these 
ranges correspond to the maximal span of gamma and 
theta/gamma oscillations for pyramidal cells. The initial 
size of the mesh has been taken as 0.01x0.01, and has 
been refined whenever necessary in order to resolve the 
fine details. It turned out that the choice of eventual 
attractor depends mostly on initial values of membrane 
potential and membrane variable of the slow cells. 

To illustrate graphically the shape of the boundary be- 
tween the attraction basins in the 12-dimensional space 
of initial values, we choose two characteristic intersec- 
tions of this boundary with coordinate planes. The first 
one (left panel of Figj9| is the plane upon which all coor- 
dinates except v p (voltage variable of the pyramidal cell) 
and ul 1 (gating variable of the left slow cell) are set to 
zero. The plot shows two distinct regions of theta oscil- 
lations separated by a strip which belongs to the basin of 
the theta/gamma state. Both borderlines display rather 
mild variations. Notably, the initial value of ul 1 should 
be sufficiently large in order to excite the theta oscilla- 
tion: initial strong contrast between the OLM-cells (re- 
call that ul 2 =0 upon this plane) would facilitate the on- 
set of the large phase difference between them, imminent 
for the theta state. 

The second coordinate plane corresponds to the case 
when only the voltage variables of the OLM-cells are ini- 
tially excited whereas the rest of the variables starts from 
zero values. In this projection (right panel of FigEJ) , the 
shape of attraction basins is much more intricate, with 
characteristic cusps, islands and narrow bridges between 
them. Here, again, we see that proximity between the 
initial states of the slow cells results, as a rule, in the 
eventual onset of the theta/gamma rhythmic pattern. 



Experimentally relevant is the location of the border 
between the basins in terms of the membrane voltage of 
slow cells. Numerics shows that this border (not pre- 
sented graphically) looks quite simple: within the stud- 
ied range of v p there exists a threshold value of volt- 
age vl,! 2 ~ 0.14 below which the unit oscillates in the 
theta/gamma rhythm whereas above it the theta rhythm 
has been registered. 




-2-1 1 2-2-1 1 2 

V P "L, 

FIG. 9: (Color Online). Projections of attraction basins of 
coexisting rhythmic patterns. Gbl x = 0.06, Gbl 2 = 0.03, 
Glp — 0.8. Dark grey: attraction basin of the theta/gamma 
rhythm. Light grey: attraction basin of theta-oscillations. 
(a): initial values for voltage of the pyramidal cell vs. mem- 
brane variable of the left slow cell; the rest of variables is set 
to zero, (b): initial values for membrane variables of two slow 
cells; the rest of variables is set to zero. 



F. Variation of asymmetry between BL-synaptic 
connections 

In order to elucidate the possible role of asymmetry in 
the pattern of synaptic connections, we performed calcu- 
lations at different values of the ratio Gbli/Gpl 2 . De- 
creasing this ratio from 2 to 1, we have detected quali- 
tative distinction only at very low values of conductance 
Glp; otherwise, the results were largely similar to the 
ones described in the previous subsections. 

In the symmetric case Gbl^_ = Gbl 2 i the system pos- 
sesses the invariant "diagonal" subspace in which the val- 
ues of variables for the two slow cells coincide (hence, 
the phase shift between them is absent). Surprisingly, 
the states from this subspace are stable not only at mod- 
erate values of Glp which correspond to the attractor 
of the theta/gamma type, but at low values of this pa- 
rameter as well, where the temporal pattern belongs to 
the gamma-type. In contrast to the previously described 
gamma-oscillations, in this state the phase shift between 
the slow cells is absent; they oscillate in unison. Below 
we call this regime a "symmetric gamma-state". If the 
conductance Glp is set below 0.0016, symmetric gamma- 
oscillations are chaotic. Above this value of Glp the 
periodic gamma-state is observed; in a narrow interval 
close to Glp=0.036 its regular pattern is transformed 
(via the mixed-mode) into the periodic theta/gamma os- 
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cillation. Outside the diagonal subspace there exists an- 
other, "asymmetric" oscillatory state which has a shape 
of gamma-oscillations with a phase difference of 27r/3 be- 
tween the slow cells and is akin to the gamma-state de- 
scribed in the Sect. IIII Al 

Accordingly, at sufficiently low values of Glp, two 
qualitatively different types of gamma-oscillations coex- 
ist: a symmetric one and an asymmetric one in which 
there is a phase lag between the slow cells. The asym- 
metric type, in its turn, consists of two different limit 
cycles: the state in which the phase of the cell L\ is 27r/3 
ahead of L2, as well as its mirror counterpart in which 
L 2 is "leading". 

Depending on the initial conditions, the symmetri- 
cal network can display either of these three gamma- 
patterns. Being structurally stable, all three limit cy- 
cles survive introduction of weak asymmetry between the 
connections; of course, the former symmetric state devel- 
ops in this case a small phase shift between the slow 
cells. Increase of Gbl ± at constant value of Gbl 2 ini- 
tially results in the growth of the attraction basin of the 
oscillation with leading L\ at the cost of the basins of 
two other gamma-states. Further growth of GbLx brings 
about saddle-node bifurcations in which at first the for- 
mer symmetric gamma-state and then the gamma-state 
with leading L2 disappear. 

For example, at Gb£ 2 =Gz,p=0.03 all three gamma 
patterns coexist in the range of values of Gbl x between 
0.03 and 0.03774. The saddle-node bifurcation of peri- 
odic orbits on the right border of this interval destroys 
the former symmetric pattern. Of the two remaining 
states with phase shift close to 2tt/3, the pattern with 
leading L 2 disappears in the saddle-node bifurcation at 
Gbl x =0.04732: beyond this value of Gbl 1 only one pat- 
tern of the gamma-type can be observed. 

Let us return to characterization of the symmetric net- 
work with Gbl x = Gbl 2 =0.03. Except for the "unison" 
symmetric gamma-state, the appearance and properties 
of the gamma- and theta-oscillations display only small 
quantitative difference from the case Gbl x 7^ Gbl 2 de- 
scribed in the preceding sections. Two paragraphs above 
we have characterized solutions which at low values of 
Glp look like a gamma-oscillation with a phase shift be- 
tween the slow cells close to 27r/3. At around Glp=0.073 
this pattern is rapidly transformed (via the mixed-mode 
oscillation) into the theta-rhythm in which the slow cells 
oscillate in the antiphase; stable theta-oscillations persist 
in the whole studied range Glp < 3. Similarly to the 
asymmetric case, at sufficiently high values of conduc- 
tance Glp, theta/gamma-state has not been observed, 
and the theta-oscillations are the only attractor of the 
system (cf. Fig. [6]). In absence of symmetry between 
BL12, the periodic theta/gamma state is eliminated in 
the course of the saddle-node bifurcation. In the sym- 
metric case, this bifurcation is replaced by the inverse 
pitchfork bifurcation; for Gbl x = Gbl 2 =0.03 this event 
takes place at Glp=2.0I4. As a result, the periodic so- 
lution of the theta/gamma type persists at higher values 



of Glp as well, but is attracting only in the invariant 
"diagonal" subspace ul 1 = ul 2 and vl 1 — vl 2 - A weak 
violation of the symmetry in initial conditions leads to 
creation and growth of the phase shift between the slow 
cells, which eventually destroys the theta/gamma oscil- 
lations and replaces them by the theta-rhythm. 

Concerning the border between attraction basins in 
case of hysteresis between theta and theta/gamma 
rhythms, merely the quantitative shift due to variation of 
the ratio Gbl x /Gbl 2 has been observed; the shape of the 
basins has been largely preserved. Compared to Fig.[9la) , 
in the symmetric case GbLi = Gbl 2 the attraction basin 
of the theta/gamma state is slightly broader. In contrast, 
projection of attraction basins for membrane variables of 
two slow cells (analog of Fig. HJb), which is symmetric 
with respect to the diagonal) displays a slight narrowing 
of the theta/gamma region. Besides, we have observed 
that symmetry in the conductances GbLx 2 lowers the 
threshold values of the membrane potentials vl 12 , which 
are necessary for the onset of the theta rhythm. 



IV. DISCUSSION 

In the preceding sections we have presented a minimal- 
istic model for a working unit of neuronal network in the 
hippocampus. Below we briefly discuss a few important 
aspects of this model. 



A. Minimality of network 

Arguably, the size of the ensemble cannot be further re- 
duced without introducing drastic distortions into collec- 
tive dynamics. For example, omittance of one of the slow 
L cells would disable the theta rhythm. As seen in Fig|5J 
in this regime two L cells are oscillating in the antiphase, 
and the fast pyramidal cell exhibits two spikes within a 
complete oscillation cycle: one spike comes shortly be- 
fore the maximum of L\ and the other one precedes the 
maximum of Li- As soon as one slow partner is removed, 
this rhythmic pattern becomes impossible. 

Further, our numerical experiments show that if the 
basket cell is removed from the configuration, the rhyth- 
mic patterns persist but their dependence on the strength 
of the coupling between the remaining pyramidal cell and 
the slow cells is weakened. If the configuration of con- 
nections is symmetric, the system is restricted to only 
one scenario of the onset of rhythms which is mostly de- 
termined by the choice of initial conditions. This seems 
to contradict to the experimental evidence which con- 
firms that switching between the regimes depends on the 
coupling strength Q- In the attempts to counteract the 
removal of the basket cell by asymmetry in the connec- 
tions from the pyramidal cell to OLM-cells we observed 
that the region of existence of the theta/gamma regime 
is significantly reduced. 
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Thus, we suppose that presence of the basket cell al- 
lows the neuronal network not only to "maneuver" be- 
tween different scenarios of emergence of all rhythms 
but also to expand the region of existence of the 
theta/gamma regime. 

The role of asymmetry in this context is twofold. On 
the one hand, as shown in Sect. IIIIF1 the symmetric case 
possesses more attractors, which enhances the flexibility 
of the network. On the other hand, reduction of variabil- 
ity in the strongly asymmetric network means increase 
in the robustness, which can also be useful in certain sit- 
uations. 



B. Possible implications for larger networks 

We have started this section by defending our choice 
of the minimal model. Here we turn to the "inverse" 
question: are patterns, generated by this minimal unit, 
recognizable on the background of large hyppocampal 
networks? Numerous experiments in the hyppocampal 
areas CA1 and CA3 unambiguously confirm existence of 
all described oscillatory states: the gamma, the theta 
and the theta/gamma [1, HI, 0|- Furthermore, in- vitro 
technique of directional cross-sections in the hyppocam- 
pus has allowed to single out the patterns: in longitudinal 
sections mostly the theta rhythm is recovered, whereas 
oscillations of gamma and theta/gamma are typical for 
transverse and "medium" (so-called coronal) slices. This 
has to do with morphology of the participating cells. 
The OLM cells are arborized (possess dendritic struc- 
ture) mostly in the longitudinal direction. In longitudinal 
slices synaptic connections from these cells to the pyra- 
midal and basket ones are at the strongest; this results in 
the synchronous theta rhythm in the area CA3 [f| 0, H(| ■ 
In other directions arborization is much less pronounced; 
hence, the relative role of the OLM cells is weaker, and 
one observes mostly the gamma and theta/gamma oscil- 
lations produced by interaction of pyramidal and basket 
cells. 

Other kinds of cells contribute to shaping and mainte- 
nance of rhythmic patterns as well: the study [To| lists 
at least 21 classes of inhibitory interneurons which par- 
ticipate in the oscillations. Anyone of them could be a 
part of the minimal module. However, high synchrony 
between inhibitory interneurons responsible for gamma- 
rhythm allows to restrict modeling to just one type: 
we have chosen the basket cells, whose interactions with 
pyramidal cells are well investigated experimentally. On 
adding the OLM cells which are required for the gener- 
ation of the theta state, we get by with just three types 
of cells in a module. 

How do such "modules" interact with each other in 
a large-scale network? What cells mediate the inter- 
actions? Which type of connection - excitatory or in- 
hibitory - should be considered? From the point of view 
of morphological properties, briefly sketched above, it 
seems plausible that modules communicate with each 



other through OLM cells. The pyramidal and basket 
cells have been shown to be involved mostly in locally 
synchronized gamma oscillations Q . Within this conjec- 
ture, OLM cells of a module can inhibit pyramidal cells 
in the neighboring modules 30]. As for interactions be- 
tween OLM and basket cells in large networks, reliable 
experimental data are still missing. 

Spatial-temporal patterns in big networks depend on 
the strength of connection between modules. For the 
case of strong coupling there is experimental evidence 
of phase theta waves [31| . In contrast, moderate cou- 
pling leads to spatial-temporal gamma oscillations which 
modulate slow theta patterns (clusters) [331 ■ Notably, 
such patterns are sensitive to the phase shift between 
the OLM cells. Finally, if the interaction is sufficiently 
weak, the modules oscillate independently of each other. 
In this case, rhythms found in a single module, persist, 
and presence of neighbors results, mostly, in shifting the 
transitions between them. 



C. Multistability 

Hystcretic properties, akin to the described effect of 
multistability, have been registered in experiments with 
subthalamic neurons [l3| . In these experiments applica- 
tion of different bias currents and shifting of the baseline 
membrane potential resulted in the onset of such different 
regimes as tonic firing, rhythmic bursts or silent upstate. 
According to the experimentalists, multistability can be 
controlled by dynamics of ionic channels [13j . 

In our model, depolarization of initial membrane po- 
tential of L-cells results in the onset of the theta/gamma 
regime, but stronger depolarization switches the latter to 
pure theta oscillations. 

To our knowledge, multistability of the discussed kind 
has not yet been registered in CA3, either in experiments 
or in existing theories. (Metastable states in large sec- 
tions of CA3, experimentally found in [l4j], refer to al- 
ternating activity of different groups of neurons, and not 
to different rhythmic patterns). In particular, the au- 
thors of the more detailed model described in [3] do not 
report on coexistence of attractors. Transitions between 
rhythms can be related to physiological state of a cell. 
The OLM cells are known to possess two specific ionic 
currents responsible for the onset of theta rhythm 0, Ht| . 
One of these currents is activated by the strong hyper- 
polarization, whereas the latter, in turn, is switched on 
by depolarization. These states correspond, respectively, 
to the top and bottom stripes in the left panel of Fig. |H1 
therefore in the corresponding ranges of potentials theta 
regime is established. For the initial conditions from the 
"middle" stripe these currents are modest, hence theta 
rhythm cannot develop, and theta/gamma rhythm is ob- 
served instead. This interpretation requires, of course, 
an experimental testing. 
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